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Abstract 

Density perturbations generated from inflation almost always have a spectral index n s which runs 
(varies with the wavelength). We explore a running spectral index scenario in which the scalar spectral 
index runs from blue (n s > 1) on large length scales to red (n s < 1) on short length scales. Specifically, 
we look for a correlation between the length scale at which n s — 1 = and the length scale at which 
tensor to scalar ratio T/S reaches a minimum for single field slow roll inflationary models. By 
computing the distribution of length scale ratios, we conclude that there indeed is a new approximate 
consistency condition that is characteristic of running spectral index scenarios that run from blue to 
red. Specifically, with strong running, we expect 96% of the slow roll models to have the two length 
scales to be within a factor of 2, with the length scale at which the tensor to scalar ratio reaching a 
minimum longer than the wavelength at which n s — 1 = 0. 



1 



I. INTRODUCTION. 



It is currently widely accepted that inflationary cosmological scenarios offer most promising 
explanations for the initial conditions for structure formation in our universe. Almost all infla- 
tionary scenarios predict that the primordial density perturbation spectrum deviates slightly 
from a power law and is dominated by scalar density fluctuations. |33( Typically, the scalar 
density perturbation spectrum is parameterized as 

p.s(k) = Au^r^- 1 (i) 

where n s (k) m 1 is a scale dependent function which is usually called the running (scalar) 
spectral index (running refers to the change in the spectral index as a function of wavenumber 
k). Combining observations of CMB, galaxy surveys, and Lya forest, there have been claims 
for evidence of strongly running spectral index (see e.g. |2j), but at the moment, combined 
data set favors no spectral index running (e.g. |3j). Nonetheless, a significant running of the 
spectral index is still a debatable possibility that will be settled by future experiments. 
According to the CMB and 2dFGRS galaxy survey data jj, [| , 

-^-(k = O.OSMpc' 1 ) = -0.03l8£g • (2) 

Furthermore, as pointed out by 0, within the context of single field inflationary models, there 
is some indication that the spectral index quantity n s — 1 runs from positive values (blue) on 
long length scales to negative values (red) on short length scales (positive to negative within 
about 5 e- folds). 

It is well known [6] that one robust check of slow roll inflationary scenario is what is usually 
referred to as the single field self-consistency condition 

n T (k) = -2e (3) 

where is the spectral index of tensor perturbation power spectrum parameterized as 
Pt oc k UT and e is the slow roll parameter characterizing the tensor to scalar power spec- 
trum ratio. In 0, it was pointed out that if running occurs from blue to red, then there may 
be another approximate consistency condition that we may observationally aim at checking 
regarding inflationary scenarios. Namely, within the single field slow roll scenario, there should 
be an approximate coincidence between the length scale k\ at which 

n s {h) -1 = (4) 



and the length scale k 2 at which the tensor to scalar ratio reaches a minimum, i.e. 

e'(k 2 ) = 0. (5) 

If true, this kind of new consistency condition would be important because currently there is 
only a very limited number of observational consistency checks that could support the picture 
of inflationary origin of density perturbations. 
The degree to which 

h « k 2 (6) 

depends on a) how fast n s (k) runs and b) the slow roll parameter e (this is true at least in the 
leading slow roll order approximation). Hence, it is not clear how compelling the expectation 
of this coincidence is. For example, does Eq. (0) occur only for 1 out of 10,000 slow roll 
models of inflation with running spectral index or does it occur for 9000 out of 10,000 slow 
roll models of inflation? Furthermore, because of the strongly running behavior of the spectral 
index, it was not clear that the leading order presented in [?| was robust. Hence, in this 
paper, we quantify the degree to which one would find the approximate consistency condition 
a compelling consistency check of single field inflationary scenario. Furthermore, we check that 
the qualitative expectation based on first order slow roll expansion is also valid at second order 
slow roll expansion. 

In the spirit of 0, we examine a grid of slow roll models and ask how the distribution of 
ki/k 2 changes as a function of increasing the requirement of running spectral index. Expressed 
in terms of difference in the e- folds of inflation (expressed as k\/k 2 = exp(AiV)), we find AN 
has a mean of about —0.34 with a 2cr width of 0.5 when a strong running condition is imposed 
as when both 0.15 < n(k = 0.002MPC" 1 ) - 1 < 0.3 and -0.3 < n(k = IMpc" 1 ) < -0.15 are 
satisfied. Hence, slow roll inflationary scenarios predict that there should be a coincidence of 
k\ and k 2 within a factor of two if the spectrum runs from blue to red, and the mismatch is 
such that k\ > k 2 . Furthermore, since the 2cr width of the distribution without strong running 
condition is about 1.2, we find quantitative evidence that if the spectral index runs strongly, 
there is an increase in the coincidence as expected. 

There are two main caveats to our results. As discussed in Section ITVl because the B-mode 
polarization sensitive to tensor perturbations is expected to peak on relatively large angular 
scales where it is measurable with current forecasts P, Q] while the minimum of tensor to scalar 
ratio lies on shorter angular scales, measurements testing the proposed approximate consistency 
condition may be difficult. If strong running is relevant to cosmology, ingenuity of scientists 



in the future will hopefully overcome this obstacle. The second is that the grid of models that 
we choose are sampled with equal weight for the initial conditions of the slow roll equation. 
Although this is commonly found in the literature say , and does not seem to be an 
unreasonable way to sample the possible set of models, this kind of probability measure does 
not have any justification from first principles. A similar ambiguity of measure plagues most 
"landscape" or anthropic principle arguments (e.g. [l| ). 

This paper is organized as follows. In the next section, we analytically compute the coin- 
cidence in terms of slow roll parameters. Following that we present numerical results. Mea- 
surement prospects are then discussed in Section ITVl In the appendix, we collect some useful 
slow roll formulas used in the derivations in the paper. Throughout this paper, we use the 
convention M p = 1/ y/G^. 

II. ANALYTICAL APPROXIMATION 

Let us see why there should be an approximate consistency condition as described in the 
introduction. In |2f , it was shown that within the context of a single real scalar field <fi slow roll 
inflationary model with a potential V(0), one can write 

where ey = ^(^V^) 2 i s the usual inflationary slow roll parameter in terms of the potential, 
the upper (lower) sign is for V'((p) > (V'(4>) < 0). Since the tensor to scalar ratio is given by 
Pg/Pn — 16ey, the vanishing of the right hand side of Eq. corresponds to the scale at which 
the tensor to scalar power reaches an extremum. Now, the observation that was made was that 
if the spectrum runs from blue to red, then ±e' v will be positive initially, but as n s — 1 becomes 
negative to the point of canceling the 2ey, e' v will vanish. As long as ey is small and/or n s — 1 
runs strongly negative, the scale at which n s — 1 crosses zero will coincide with the extremum 
of ey which in turn corresponds to the extremum of the tensor to scalar ratio. 

In this section we will extend the previous derivation to second order in slow roll and confirm 
the qualitative validity of the previous results. This is reassuring given that the validity of 
leading order results was questionable in light of strongly running behavior. To accomplish 
this task, we will work with Hubble flow slow roll expansion techniques (see, e.g. pj). The 
advantage of the Hubble flow approach is that it is simpler to go to higher orders in derivative 
expansion. 



Using the conventions of |lJ], we use the definitions for the slow roll parameters in terms of 
the Hubble function H (0) as follows: 
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Here, the Hubble function is defined by the Hamilton- Jacobi equations 
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which has the advantage of having a readily solvable model of power law inflation. As is well 
known, Hubble flow slow roll expansion that we summarize below is a derivative expansion in 
H((j>) about the power law inflationary model. 

We used iV as the measure of time during inflation, the number of e-folds before the end of 
inflation, which increases as one goes backward in time : 
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These equations implies a relationship of e to H and iV: 
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The evolution of the higher order parameters during inflation is determined by a set of "flow' : 
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As one can see, in general any slow roll parameter could be expressed in terms of one single 
slow roll and its derivatives. Some of these relationships are explicitly given in the appendix. 
Note that e has the interpretation through the fact that tensor to scalar power spectrum is 
given by 

= 166 (16) 

where P g and are the usual power spectrum for tensor and scalar gauge invariant metric 
perturbations. 

Let us now solve for the difference between the length scale at which e reaches an extremum 
(l/k 2 ) and the length scale at which n s — 1 crosses (1/fci). Instead of giving the length scale 
in terms of the inverse wave vector, we will express it terms of the difference in the horizon exit 
e-folds AA. To translate that into ratio of lengths, one simply has 

In A = AN = N l -N 2 + 0(\n(H 2 /Hi)) (17) 

where H 2 is the expansion rate when k 2 left the horizon and Hi is the expansion rate when ki 
left the horizon. Since \n(H 2 /Hi) is of order of the slow roll parameter during the early stages 
of inflation, one can disregard this term when interpreting AiV = Ni — N 2 as long as it is larger 
than about 0.1. Defining iVo as the point where the extremum of e occurs, i.e. ^\n = 0> we 
can linearize n s — 1 around it to solve approximately for AiV. Using the formulas given in the 
appendix, we can write 

n s (N + AN) - 1 ~ n s (A ) - 1 + ^k,AA (18) 



dn s 2e e — e — ee 



dN e 
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(19) 



n s -l = -Ar{(-3 + C)i 2 + he + 2e 2 (C + 3)1 e + e(3 - C) - 8e 3 (e + 1)}. (20) 
Solving n s (N + AN) - 1 = for AiV, we obtain 

a at C-3 2e(l + e) 

AN ana ly tica l RS - h d . ■ (21) 

The constant C = 4(ln2 + 7) - 5 = 0.0814514, where 7 ~ 0.577 is Euler's constant. As it can 
be seen from this formula AA^ clusters around the negative value —0.67 and it can be positive 
or negative. We can also observe that stronger running (associated with a higher slope for n s ) 
drifts AA toward more negative values, and that this effect is modulated by the value of e at 
A . 
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FIG. 1: From the graph of e , Pji (scalar power spectrum), and n s — 1 as functions of the number of 
e-folds N, it can be seen that the extrema of e and Pji are closed to the point where n s — 1 is zero. The 
dot dashed curve is Pji/Pji(0), the solid curve is 50e, and the dotted curve is n s — 1. The horizontal 
axis is oriented such that inflation ends at N < —3.5 (i.e. time flows right to left). 

We will see that numerical results agree well with the analytical approximation. Before 
going onto numerical analysis, we would like to note that beyond Eq. (J7|), there may be many 
more "approximate consistency" conditions associated with strong running since we are merely 
labeling various features of the slow roll equations. For example, as it can be inferred already 
from the first order formula = the power spectrum will have an extremum roughly 

coinciding with the extrememum of e . To check this numerically, we can examine an example 
of an inflationary model (as explained in the next section). As it can be seen in Fig. ^ the 
scalar power spectrum has an extremum very closed to the point where n s — 1 is zero. We could 
proceed similarly to the way we have done for Eq. ©, comparing analytical approximation 
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to numerical data for this "approximate consistency" condition, which in this case would not 
depend on the strength of the running of the spectral index, but could be another test for more 
general inflation models. 

III. NUMERICAL APPROACH 

In this section, we study our proposed approximate consistency condition question numer- 
ically. One scientific question we would like to address is, "What is the expected coincidence 
of the value k% satisfying n s (kx) — 1 = with the value k 2 satisfying e'(0fc 2 ) = 0?" The sec- 
ond question is whether a stronger negative running of the spectral index makes the expected 
\k\ — k 2 \ significantly smaller. 

To this end, we examine a particular grid (to be specified below) of slow roll inflationary 
models. In this grid of models, we make a set of "cuts" which select out subsets of the models 
with scalar spectral index running from blue to red. In this way, we can make a histogram 
of models as a function AiV (defined in Eq. (|T7j) ). The peak and the width of the histogram 
quantifies the expected coincidence of k\ with k 2 , answering our first question. We then make 
the cuts more stringent such that only those slow roll models with more strongly running 
spectral indices are included in the histogram. The degree to which the width of the histogram 
changes as we make the required running stronger answers our second question. 

Let us now explain the details which closely follows the methods of js}]. Assigning a set of 
initial conditions to Eqs. (fTo^l and integrating them specifies a flow of the slow roll function. 
Each of the initial conditions corresponds to a particular choice of inflaton potential and the 
inflaton initial conditions. We choose initial values for the parameters at random from the 
following ranges, assuming a uniform probability distribution : 

[40, 70] 
6 10~ 76 ,0.6" 
[-0.5,0.5] 
[-0.05,0.05] 
[-0.025,0.025], 

A1+1 A H = 0. (22) 
We truncated the expansion to order 5 by setting 6 Ah = 0. We calculate the values of the 
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tensor/scalar ratio r, the spectral index n, and the "running" of the spectral index dn/dlnk 
according to: [l^, [3] 

r = 10e[l-C(«r + 2e)], (23) 

n s - 1 = a - (5 - 3(7) e 2 - ~ (3 - 5(7) (re + i (3 - C) ( 2 A H ) (24) 

Note that the variable r here is defined with a different normalization than in where r « 16e 

to leading slow-roll order. A comoving scale k crossed the horizon a number of e-folds N(k) 
before the end of inflation : |2C 



^) = 62-ln^-lni^+ln^-IlnV ( 25 ) 

Here 14 is the potential at horizon exit, V e is the potential at the end of inflation, and prh is 
the energy density after reheating. Note that since slow roll inflation evolves toward decreasing 
potential, we can write 

V e = f£V k (26) 

where /& < 1 is a function of k that varies as to keep V e independent of k. Furthermore, since 
the reheating energy density must be smaller than V e , we can write 

Prh = 7 4 K (27) 

where the proportionality constant satisfies 7 < 1. Using Eq. ([25)1 . derivatives with respect to 
wavenumber k can be expressed in terms of derivatives with respect to iV as [2j| 

5iv = - (1 - £) dhS' (28) 

Hence, the running of the spectral index can be expressed to third order in slow roll parameters 
as 
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+ (-10 C + 6) ( 2 A H )e + (C - 3) ( 2 A H )a + +20 e a - 8 ( 2 A H ) + (2 C - 6) ( 3 A H )](29) 

We will now use superscripts a and b to denote quantities evaluated at or corresponding 
to the length scale k a = 0.002 Mpc -1 or k b = 1 Mpc -1 . The algorithm for generating the 
histogram of Fig. 1 can then be described as follows: 

1. A set of random initial conditions was generated for the value of the slow roll parameters 
at N a . The following constraints are imposed: 

r a < 0.5 (30) 



n? < < - 1< < 



< 



dink 

where n„ is fixed to 0.3 and nf varies from to 0.225 for different cuts with a 0.075 
increment. If these constraints are respected, proceed to step 2 otherwise go back to step 
1. 



2. Integrate the flow equation (using LSODA from ODEPACK [15]) to N and then check 
the following constraints: 

n\ < nl - 1< n\ 

where n\ is fixed to -0.3 and n h u varies from to -0.225 for different cuts with a 0.075 
increment. If these are respected go to step 3 otherwise go to step 1. 

3. Evolve forward in time (dN < 0) until inflation ends (e > 1) .If the total number of e-folds 
N from the beginning to the end of inflation is in the range [40,70] add this model to the 
ensemble of acceptable models. If not go back to step 1. 

4. Repeat step 1 through 3 until the desired number of acceptable models have been found. 

Once the set of acceptable models has been generated we solve numerically the following 
two equations: 

3*l» = ° (31) 

n s \ Nl -l = (32) 

It is important to observe that the presence of the square root of e in the definition of the 
differential operator (Eq. (fT2"j) ) causes numerical solutions to Eq. (|3*T|) to be difficult for 
small values of e, since small values of the derivative may be associated to the smallness of e 
more than to the actual presence of an extremum. In other words, even though the derivatives 
are taken with respect to N, for inflaton potentials with very small e, the derivatives become 
very small, requiring high numerical precision. For this reason and to discard inflection points, 
extrema are identified numerically using the second derivative information as well. 

Before making the histogram, we must further check that the models in the histogram are 
reasonable from an inflationary phenomenology point of view. This is done in lieu of executing 
precision fits to all available data since it has been well demonstrated by previous analyses (see 
for example |8j) that the current CMB data are not extremely constraining in terms of the 
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details of the slow roll inflationary models. Hence, our histogram should not be sensitive to 
the lack of precision in our cuts. On the other hand, given that the combined fits including 
SDSS data as carried out by |3J] prefer small running, the strongly running cases of the current 
analysis may have been statistically disfavored had precision fits to combined data been made. 
Given that most of the numerically generated models do not have strongly running spectral 
index for the weakest cut, the results of the weakest cut should be robust with respect to 
imposing additional fit constraints, while the results of the stronger cuts should be interpreted 
with appropriate caution. 

Let us now write explicitly the constraint equations. Since the WMAP analysis [3] gives 

p((k a = 0.002MPC" 1 ) = (2.95 x i(T 9 )(0.75 ± 0.09) , (33) 

we use the approximate formula for the power spectrum 

C _ 8V k 

^ " 3eA# ( j 

to impose a constraint on e and V k . The requirements of having sufficient energy density at the 
end of inflation to reheat the universe to a temperature consistent with big bang nucleosynthesis 
gives 

(10- 2 GeV) 4 < prh <V e < V k . (35) 

We also know from observations that r a < 0(1) and that Eq. (|33J) is true. Substituting V k 
obtained from Eqs. (|25p. ([26)1. and (}2~Tj) into Eq. and solving for e we obtain 

6 x 10~ 76 < e a < 0.05 (36) 
e a = e^~^y(ka) (37) 

where y(k a ) = f ka r )~ A ^ and the upper bound comes from Eq. (|3T)|) and r = lOe (again note 
that the normalization used in this paper is given by Eq. (|2*3*j) . and specifically is not r = 16e). 
Consider how Eq. ()37|) gives a constraint. Note that for any given flow trajectory, because of 
Eqs. (|3^j) and (jHIj), V ka is fixed. Furthermore, the flow trajectory itself determines f ka (the 
value of the potential at the end of inflation). Hence, the only adjustable parameter is 7 which 
determines the reheating temperature. For every flow trajectory with N(k a ) fixed, Eq. (J37)) 
then provides constraints on the initial conditions of inflation. For example, consider a flow 
trajectory having N(k a ) = 70. Eq. (|3*7j) gives 

ft a < 0.05exp(-4[iV(A; a ) - 61.8]) ~ 10" 16 (38) 
11 
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TABLE I: Basic distributional characterization of Fig.|2J m<i = \J ((AN) 2 ) characterizes the deviation 
from 0. P((AN) ± 2<taa0 is the fraction of the models within 2gan of (AN). N m d is the number of 
models analyzed for the given range of the n s . 



This is a significant constraint since this says that the energy density at the end of inflation is 
at least 10~ 16 smaller than the the energy density at 70 e-foldings before the end of inflation. 
The physical reason for Eq. (|36|) is that upper and lower bounds on the potential (coming from 
tensor perturbations and reheating temperature, respectively) directly translates into upper 
and lower bounds e because of Eqs. and (|54~j) . The physical origin of Eq. (|37j) is that 
any mismatch in the inflationary stretching of the wavelength and physical k/a = 0.002Mpc~ 1 
must be compensated by postinflationary expansion which is related to the potential energy 
at the end of inflation which in turn is related to e through Eqs. (jH4j) . and 1)26)1 . In 

practice, Eq. ()38j) does not play a significant role in the numerical exploration since the number 
of e-foldings is generically small when there is significant running. 

Our final histogram can be seen in Fig. Eland its distributional characterization can be seen 
in Table |l] To have a measure of how AiV clusters around zero we followed this procedure: a) 
we have divided AiV in the range [—1,5.25] into 25 intervals of width 0.25 each containing n; 
models and then calculated for each of these intervals the average (AiV)j for the models in that 
interval and the standard deviation <j ANi = ; b)we calculated m 2 = S?£ x ^a^* , a weighted 
second order moment around of the 25 averages of (AN)i for the intervals denned above, 
using the inverse of the corresponding variance a ANi as weight. We also similarly calculated 
the weighted standard deviation a m2 . The total number of slow roll models we considered to 
construct the histogram was around 10 7 , and the number of models fulfilling the cuts is shown 
in the N m d column of Tbl. HJ 

Notice in Table Ulthat both |(AiV) ±<taat| < 1 and |m 2 ±cr m2 | < 1 regardless of the strength 
of running. The column labeled " P((AN) ± 2a an)" quantifies (at least with respect to our 
choice of model sampling) the extent to which we should expect a coincidence between k\ and 
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FIG. 2: Frequency histogram (hist.) for the difference between solutions for different running of the 
spectral index, i.e. different n a (spectral index at wave vector k a = 0.002Mpc~ 1 ) and n& (spectral 
index at wave vector k b = IMpc -1 ). Thick solid hist, gives the results for { < n a — 1 < 0.3 , 
-0.3 < n b - 1 < 0}; dashed hist, for { 0.075 < n a - 1 < 0.3, -0.3 < n b - 1 < -0.075}; the green 
hist, for {0.15 < n a - 1 < 0.3 , -0.3 < n b - 1 < -0.15}; the solid hist, for {0.225 < n a - 1 < 0.3 , 
—0.3 < n b — 1 < —0.225}. The area under each histogram is normalized to unit area. 



k2- i.e. we should expect with roughly "95% confidence" that the length scale for which the 
tensor to scalar ratio reaches a minimum will coincide with the length scale for which the 
n s — 1 = up to a factor of about 2 and k\ > k% (i.e. wavelength at which the tensor to 
scalar ratio reaches a minimum is longer than the wavelength at which n s — 1 goes through a 
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zero). As far as the effect of different strengths of running is concerned, one sees in Table H] 
that stronger running corresponds to a larger magnitude of | (AiV) | and a more negative value 
of (AiV). Although "statistically" marginal, the smaller catv means that strong running seems 
to narrow the distribution of AiV. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the robustness of the approximate condition suggested by [3| • 
We have quantified distributionally the extent to which one should consider the approximate 
consistency condition to be a prediction of single field slow roll inflation. We find that if the 
spectrum runs from blue to red, we should look for an approximate coincidence between the 
wave vector k\ satisfying n s {k\) — 1 = and the wave vector k 2 satisfying e'fa) = 0, up to a 
factor of about 2. According to a WMAP analysis Q|, this should occur at k/ao ~ 0.02Mpc -1 . 

Given that the consistency condition is only approximate, even if one finds that observations 
contradict the consistency condition (for example, if observations deduce k\jk 2 < 1) it will be 
very difficult to make any rigorous conclusions about the models of inflation. Furthermore, 
given that we have sampled each of the initial condition ranges with equal weight, it is not clear 
what the distribution that we computed has to do with the real world. (A similar argument 
can be made about much of the current literature that use similar sampling 2, |^, 11, 12 1.) 
Nonetheless, if this set of initial sampling turns out to be a good approximation to the real 
inflationary model distributions coming from a more fundamental theory, tests coming from 
this kind of approximate consistency condition may be useful. To turn this around, there may 
be a way to classify the vacua of more fundamental theories according to whether they satisfy 
roughly equal probability sampling of the slow roll initial condition or whether they prefer a 
very specific distribution which makes the approximate consistency condition more (or less) 
stringent. 

Before concluding, let us briefly consider the measurement prospects for the approximate 
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consistency condition. According to current ideas of anticipated measurements 
B-mode polarization observations are necessary to reliably infer the tensor perturbation ampli- 
tudes. As is well known, B-mode polarization has the advantage over the E-mode polarization 
in that at the last scattering surface, the B-mode polarization cannot be generated by scalar 
perturbations alone because to leading order, the scalar perturbations do not generate nonzero 
{/ = 2, m = ±2} amplitudes in the temperature anisotropy while the tensor modes can. Hence, 
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in principle, assuming no large contamination from vector perturbations, B-mode polarization 
measurement offers a direct way to infer tensor perturbation amplitudes. 

The difficulty with measuring the B-mode polarization spectrum coming from gravity waves 
is that because its amplitude is proportional to the local quadrupole anisotropy induced by the 
time derivative h (where h is the amplitude of the gravitational wave) localized at the last scat- 
tering surface (due to the requirement of Thomson scattering), its amplitude is suppressed at 
least by order k/ (anHpi) on long wavelengths with respect to the temperature anisotropy (where 
Hfj and cir are the Hubble radius and the scale factor at recombination, respectively). Further- 
more, since gravitational wave mode amplitudes inside the horizon are diluted by the expansion 
of the spatial volume, the transfer function should behave approximately as t where 

ao is the scale factor today and a^/ao ~ 10 3 is the redshift to the last scattering surface. (Here, 
we have considered only the modes that entered the horizon after matter domination, i.e. 



k/ao < Hq^I ao/ a eq .) Hence, the B-polarization temperature perturbation amplitude can be 
written as 



k^Q B ,(k) ~ 10- 1 k R 1 Jl^l ) x U^k/(a H )) (39) 
where is the power spectrum of the tensor perturbations. (For discussions of analytic treat- 



ments of CMB polarization, see e.g. 



27|,|2i 



29j, |30j|.) This peaks at k/a^ = ifoy a o/ a * ~ 
l/(100Mpc), and even there, the amplitude is suppressed by about 10~ 2 relative to the tempera- 
ture anisotropies if the tensor perturbation amplitude is about the same as scalar perturbations 
amplitudes. 

To make the situation slightly worse, k% is expected to be at around k\/ao ~ l/(50Mpc). 
(Note that to map the wave vectors k/ao to the multipole moment number I, one can use the 
approximate formula — ~ -^rHn rs ^r-r — •) This means that to actually measure ki/ko, 

11 a 3.5 u 15,000 MpC ' J 1 ' 

one must go to short scales (I ~ 300) while the B-mode polarization measurements prefer large 

nn 

length scales (I ~ 90; for more accurate plots of examples, see e.g. |9j,[10j). Hence, experimental 
confirmation of the approximate consistency conditions will be a challenge. 

There is one more obstacle that makes it difficult for tensor perturbation amplitudes to be 
extracted from B-mode polarization measurements. This is due to the fact that E-modes can 
be converted into B-modes through gravitational lensing Hence, obtaining a tensor 

spectrum to check the approximate consistency condition requires success in extracting the 
tensor perturbation contribution even at relatively short length scales, around where the lensing 
contribution most likely dominates over the gravitational wave signal. This contamination from 
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gravitational lensing can be subtracted out, at least in principle, if the lens distribution can be 
accurately deduced. 

Despite the difficulties, we are optimistic that the ingenuity of researchers in the future will 
allow a test of the approximate consistency conditions considered here. 
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APPENDIX A: SLOW ROLL FORMULAE 



Here we report some formulas which can be used to express slow roll parameters in terms of 
e and its derivatives respect to the field or the operator Denoting ^ = e and ^ = e : 

4 e 3 / 2 7r + e i/7rm 
V = 477? 

-4 e 3 / 2 vr + e y^Mp 



a 



2 A H 



V 
a 

2 A H 



16 7re 3 + 12e 3 / 2 v^Fe' M p - e' 2 M p 2 + 2e" M p 2 e 



7T (2e 2 + e) 



2tt e 
e-2e 2 



2e 4e +3e 2 e-e 2 + ee 
2? 
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